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A primordial magnetic field (PMF) can affect the evolution of density field fluctuations in the 
00 ', early universe. In this paper we constrain the PMF amplitude Bx and power spectral index ub 

■ by comparing calculated density field fiuctuations with observational data, i.e. the number density 

' fluctuation of galaxies. We show that the observational constraints on cosmological density fluctu- 

(N ■ ations, as parameterized by erg, lead to strong constraints on the amplitude and spectral index of 

^ ! the PMF. 

^ ■ PACS numbers: 98.62.En,98.70.Vc 

m 

I ' 

O ' explanation for such magnetic fields in galactic clusters is the existence of a primordial magnetic field (PMF) of 
^ I order 1 nG whose field lines collapse as structure forms. Therefore, recently, the origin and amplification mechanism 
^ . of the PMF in the scale of galaxy cluster have been proposed and studied intensively by a number of authors 
L_j ; 0, B 0, H, a E [H [3 [H, El • The origin and detection of the PMF is, hence, a subject of considerable ^terestjn 
modern cosmology. Moreover, the PMF could influence a variety of phenomena in the early universe 
such as the cosmic microwave background (CMB) fl^, 'To', '2^, '2T, '55, 111, [11, 111 111 113, III, 
and the formation of large-scale structure (LSS) ^37., ,38.., .39,,, .40,. 41,]. 
' If dynamically significant large-scale magnetic fields were present in the early universe, they would have affected 
' the formation and evolution of the observed structure. Thus, some signatures of the existence of a PMF should be 
. apparent in the presently observed cosmic structure. 

In this regard, the alternative normalization parameter as is of particular interest. It is defined [4^ as the root-mean- 
square of the matter density fluctuations in a comoving sphere of radius 8h~^ Mpc. It is determined by a weighted 
' integral of the matter power spectrum. Observations which determine ag provide information about the physical 
processes affecting the evolution of density-field fluctuations and the formation of structure on the cosmological 
ILJ ' scales. The mechanisms by which a PMF can affect the density field fluctuations on cosmological scales has been 
• [ described in our previous work [sj. Of course, as is also affected by the presence of a PMF. In this article we show 
that by considering the effect of a PMF on as and comparing theoretically estimated values for cg with the observed 
?H range, we can obtain not only insight into the underlying physical processes of density field fluctuations in the presence 
. 5^ . of a PMF, but also place constraints on the amplitude and spectral index of the PMF. 



I. INTRODUCTION 



Magnetic fields have been observed [U, H, H, 0] in clusters of galaxies with a strength of 0.1 — 1.0 /i G. One possible 




II. THE MODEL 



We use the isocurvature magnetized initial conditions with adiabatic relations for the fluids evolution of primary 
density perturbations and in the presence of a PMF. For the present purposes we fix the cosmological parameters to 
those of the best-fit flat ACDM model as given in Ref. i.e. h = 0.719, flth'^ = 0.02273, Oc/i^ = 0.1099, ns = 0.963, 
and Tc — 0.087, where h denotes the Hubble parameter in units of 100 km s^^Mpc"^, fib and flc are the baryon and 
cold dark matter densities in units of the critical density, ns is the spectral index of the primordial scalar fluctuations, 
and Tc is the optical depth for Compton scattering. We use natural units c = h = 1. 
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A. Primordial Magnetic Field 



Before recombination, Thomson scattering between photons and electrons afong with Coufomb interactions between 
electrons and baryons were sufficiently rapid to ensure that the photon-baryon plasma behaved as a single tightly 
coupled fluid. Since the trajectories of plasma particles are bent by Lorentz forces in a magnetic field, photons are 
indirectly influenced by the magnetic field through Thomson scattering. The energy density of the magnetic field 
can be treated as a first order perturbation upon a flat Friedmann-Robertson- Walker (FRW) background metric. In 
this linear approximation [31, the magnetic field evolves as a stiff source. Therefore, we can discard all back reaction 
terms from the magnetohydrodynamic (MHD) fluid onto the field itself. 



B. Power Spectrum from the PMF 



During the epochs of interest here, the conductivity of the primordial plasma is very large, and the PMF is "frozen- 
in" to a very good approximation [23. Furthermore, we can neglect the electric field, i.e. E ^ 0, and can decouple the 
time evolution of the magnetic field from its spatial dependence, i.e. B(x, r) = B(x)/a2 for very large scales. This 
leads to the following simplified electromagnetic energy-momentum tensor, 

T [em](x,t) = ^ , (1) 
r'°[EM](x,T) -T°'=[em](x,t) =0 , (2) 
-r^EM](x,r)=< = ^{2i?'(x)i?'=(x)-<5^'=i?(x)2} . (3) 

We assume that the initial PMF is statistically homogeneous, isotropic and random. For such a magnetic field, the 
power spectrum can be taken as a power-law P{k) oc fc"^ [25l l3l| where ub is the spectral index which can be either 
negative or positive depending upon the physical processes of the field creation. From ref. [25|| , a two-point correlation 
function for the PMF can be defined by 



(n \nB+S d2 

B\k)B^'{k')) = i-2^_^-^fc"Bp^^(fc)5(fc-fc'), k<kc, (4) 



where 



P^\k) = S^'-^ ■ (5) 

Here, Bx is the comoving mean magnetic-field amplitude obtained by smoothing over a Gaussian sphere of comoving 
radius A, and kx = Itx jX (with A = 1 Mpc in this paper). Hereafter, we work in k-space and denote all quantities by 
their Fourier transform convention 

i.,.)^/.>„xp,...x)f(x). (6) 

The cutoff wave number kc in the magnetic power spectrum is defined by [45l . , 



T{p+p) Jo T' ^ < ^de 



fc^"-"^(r) = <^ 47r{p+p) Jo T' ' ^ 'dec (7) 



where is the photon mean free path, and Tjec is the conformal time at the epoch of photon-baryon decoupling. 

We obtain power spectra for the PMF energy density and the Lorentz force for the scalar mode, respectively, as 
follows 

1 



|£^[EM:S](fc,T)|25(fc-fe') = ^(T(fc,r)[EM:Sl]r*(fc',T)[EM:Sl]) , (8) 

and 



|n[EM:S](fc,T)|2<5(fc-fc') = ^;^-^((T(fc,T)[EM:Sl] -I^(fc,r)[EM:S2] 



(27r)3 

X (T*(A;',r)[EM:Si] - 7'*(fc',T)[EM:S2])) , 

(9) 
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where SI and S2 in the subscripts of the energy-momentum tensor denote the PMF energy density and pressure. An 
exphcit expression can be obtained for the ensemble averages which are used to evaluate the above spectra. In the 
case of a power law stochastic magnetic field we havefsil [13, l47l.[4^. 



(27r)"+s Bl 



5{k-p) 



{U ■ {k-U)Y 

k'^\k-k'\^ 

The two-point correlation function for the Lorentz force is given by 

1 f (27r)"+^ B'^ 1 ^ 

(T(fc,r)[EM:Sl]?'b,T)[EM:Sl]) = 23(27r)7a8 | 2/;^+^ fj^^ J 

' J dC\k- fc'|"-2 {(1 + C^)k^ - 4kk'C + 2k'^} S{k - p) 



(10) 



Here, we define C as 



C = 



= kk = 



U ■ k 

k'k 



(11) 



(12) 



Almost all previous works have set the terms which include C in the middle parenthesis to unity. In this paper, 
however, we evaluate Eq. (jlip explicitly using integration by parts. In this way we obtain the following equation. 



{T{k, T)iEM:Si]T*ik, r)[EM:Si]) 



1 J (27r 



87ra8 1 2kl''+^ r("«+^^ 



dk'k"'^+^ 
1 

k''^nB(nB + 4) 
k 



n| + Aub + 1 



kk'uBiriB + 2){nB + A) 
{|fc-A:'|"^+2 + |fc + /c'|"«+2} 



{(fc + fc')"^+^ - |fc-/c'|"«+2} 



fc'3ns(riij + 2)(nB+4) 



{{k + k') 



t\nB+2 



-\k-k 



nnB+2 



(13) 



A similar derivation leads to the power spectrum of the PMF tension as well as the power spectrum of the correlation 
between pressure and tension as follows. 



TlEM:S2]{k)T^EM:S2]W/ " 
X / dk'k''^''+^ 

3 



1 j (27r)"«+5 Bl 



{kk'){7lB+&) 

3 



{kk'fnB[nB+2){nB+A) 
{|A:-A:'|"^+*^ + (/s + fc')"^+^} 



{{k + k'Y''+'^ -\k-k'Y'''+^] 



(A;fc')2(nB + 6)(nB -f 8) 



{(fc-t-fc')"''^^ - \k-k'Y'^+'^] 



(14) 
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and 



(%M:Sl](fc)7'*EM:S2]('^)) + (%M:S2] (fe)7'*EM:Sl] ('^)) 

2 



1 



(27r 



+3 



(fcA:')2nB(nB + 2) 



{(fc + fc') 



l\nB+3 



|fc-fc'|"«+3} 



Ph'^nsinB + 2)(nB +4) 
1 

fc3fc'2nB(nB + 2)(nB+4) 
3 



k^k'^nsinB + 2){nB + 4)(nB + 6) 



{|fc-A:'|"«+'' + (fc + fc')"^+''} 

{{k + k')''^+^ -\k-k'\''''+^} 



(15) 



For this article we have constructed a numerical program, "PriME: Program for primordial Magnetic Effects" , with 
which we can evaluate the PMF source power spectrum using the numerical method described in Refs. [23l. [stI. [49| . 
Using this, we can quantitatively evaluate the time evolution of the cut off scale and thereby reliably calculate the 
effects of the PMF. 



III. EVOLUTION EQUATIONS 

We now summarize the essential evolution equations for each mode. 

For the scalar mode we obtain the following equations in fc-space [H, [13, HH US [13] : 



k^<P + 3H{^ + Hyj) = ATrGa^{E[EM:S]{k,r)-Sptot} (16) 
k\<p-4>) = -l27rGa^Z[EM:S]ik,T)-ip, + P,)a,} 

= -UnGa^ I^^E[EM:S]{k r) +n[EM:S]{k, t) - {p, + P,)a,^ , (17) 



Snp I 



^ [E[EM:S]{hT)+6HE[EM:S]{k,T)} , (18) 



= _^(l_3^,yS)_^„(S) + |P^!^_fc2^ + ^2^ 

^ ' 1 + w Spl + w 

^^,ntEM:Sl(fc,r) ^^g^ 
47rp 

where ^[em:S] = fB(P7 + P7) and w = p/p is the usual equation of state parameter. Note that for the photon 

= 40g^\ and u^^^ = k&^-^\ Massless neutrinos obey Eqs. (US]) and as written without the Thomson coupling 
term. In the continuity and Euler relations (Eqs. fTSland fTO)) for the scalar mode, we can just add the energy density 
and pressure of the PMF to the energy density and pressure of the cosmic fluids. Since the baryon fluid behaves like 
a nonrelativistic fluid during the epoch of interest, we can neglect w and SPj^^^ /6p^^\ except for the acoustic term 

Csk'^S^^^K Also, the shear stress of the baryons is negligible [ll|. Since we concentrate on scalar type perturbations in 
this paper, we do not consider the magneto-rotational instability from the shear stress of the PMF and baryon fluid 

By considering the Compton interaction between baryons and photons in equations (jlSp and (|19p we obtain the 



same form for the evolution equations of photons and baryons as in previous work [501 1 511 . 152l 153 



k\ct>-^) 



5. 



V. 



CDM 

(S) 
CDM 

(S) 

7 

(S) 
?(S) 



.(S) 



V, 



V 



AwGa^ {£;[EM:S](fe,T) - Sptot} 
-12TTGa^ {ZiEM:S]{k,T) - {p^ + P^)(J^} 



{pu + Pu)(Jv 



_1„(S) + 

3 ^ 
3 " 
4^ 



40 
40 



^7 



I (s) 



(S) 

b 



30 



a (S) 

— < 
a 



t>Ph ^ 



"b ) 



3^^2 n[EM:S](fc; t) 

4 



(20) 

(21) 

(22) 
(23) 

(24) 

(25) 

(26) 

(27) 
(28) 

(29) 



where R = (3/4)(pb//9^) is the inertial density ratio of baryons to photons, rif, is the free electron density, ut is the 
Thomson scattering cross section, and of the second term on the right hand side of equation ([26]) is the shear 
stress of the photons with the PMF. Since ne < is favored by constraints from the gravitational wave background 
[s^ . and the PMF effects are not influenced by the time evolution of the cut off scale kc for this range of tib, we can 



approximately set i?[EM:S] 



in the following analysis. 



A. Initial Conditions 



We need to specify the initial perturbations for solving the evolution equations presented in the previous section. 
We start the solution at early times when the k modes of interest are still outside the horizon, i.e. the dimensionless 
parameter fcr ^ 1. We consider only the radiation-dominated epoch since the numerical integration for all of 
the k modes of interest will start within this era. Baryons and photons are tightly coupled at this early time 
and the expansion rate is if = t^^. We derive initial conditions for all of the modes utilizing the method of 
Refs. [35l. [sol. [Sll. [S^. [ssl. [S^ . We can assume that all density fields are zero initially, since the PMF only affects the 
velocity field of ionized baryons, via the Lorentz force, and the density fields are not directly affected by the PMF. In 
the radiation dominated epoch, photons and neutrinos are important in the energy-momentum tensor. The evolution 
equations for the photons and neutrinos are 





= +40 , 


(30) 




= -^4^^+40 , 


(31) 






(32) 




= k' (^16^ - aP + fcV) , 


(33) 






(34) 
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Here, we have omitted higher muhipole moments £ > 1 for photons and £ > 2 for neutrinos. At fowest order in kr, 
the initial conditions for Eqs. (|20ll25)) are 

= - - p^^u - R.Rb + 4i^ /7^+^;f , (35) 

(S) ^ _15fl,4aB + fl.i?B 2 (37) 

4 i?^ 4i?^ + 15 ' ^ ^ 

i?^ 2i?^. 4i?^ + 15 ' ^ ^ 

V. = -20 = -2i?,— ^^^^^^ , (39) 



where 



Rv 
Rb 



Pi 



p^ 

P^ + Pu' 
£'[EM:S] 



Pi 



IV. MATTER POWER SPECTRUM 

Possible origins of the PMF have been studied by many authors, however, there is no consensus yet as to the origin 
of the PMF. Thus, we cannot know how the PMF correlates with the primordial density fluctuations. However, 
almost all previous works investigated the effects of a PMF on density perturbations under the assumption that there 
is no correlation between the PMF and the primordial density fluctuations [131 ■ However, in order to study the PMF 
effects in a more general manner, we introduce a parameter "s" which characterizes the correlation between the PMF 
and the primordial density fluctuations [35l [37l| . In the linear approximation, the power spectra of the baryon (Pb(fc)) 
and CDM (PcDM(fc)) density fluctuations in the presence of a PMF are then written, 



Pb(fc) = (^'5[b:FL](fc)'5fb:FL](fc)/ + \VpMF](fc)^*b:PMF](fc)^ 

+ 2(<5[b:FL](fc)'5[b:PMF](fc)): (40) 

PcT,M{k) = ('5[CDM:FL](fc)'5*cDM:FL](^)) + ('^[CDMiPMF] (fc)^[bDM:PMF] (fc) ) 

+ 2<^(5[cDM:FL](fc)'5fcDM:PMF](^)) ' (41) 

where we normalize the cross correlation terms with the parameter s, 

(l5[b:FL](^)<5[b:PMF](fc) ) = S\j (^[b:FL](fc)^*b:FL](^) ) ( '^[b:PMF] (fc)^*b:PMF] (fc)), (42) 
( <^[CDM:FL] (^)^*CDM:PMF] 

(fc) ^ sM 

•^[CDMiFL] (^)'5[CI3M.FL] (^) ) ( '^[CDM:PMF] (^)'5[cF)ivi,piviF] 

(fc) . 

(43) 

Here, 5a, a ^ ([b : FL], [CDM : FL]) designate the baryon and CDM density fluctuations without the PMF respectively. 
Similarly, (5,5, ^ e ([b : PMF], [CDM : PMF]) denote the baryon and CDM density fluctuations with the PMF included. 
When < s < 1, s = 0, or — 1 < s < in eqs. (j42|) and (|43p. one has positive, vanishing, or negative correlations, 
respectively. 

The square root of the power spectrum functions for the Lorentz force H[em:S](^j ''') in Eq. (j29p does not specify the 
sign. In other words, there is no information as to whether the magnetic pressure or the tension is dominant, and 
whether the directions of forces from them are the same or different. However, such information should be taken into 
account. 
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The Lorentz force term in Eq. (|29p can be divided into two terms, the magnetic pressure and the tension. Their 
ampHtudes are given by Eqs. (|13p and p^ . respectively. By comparing those equations, one can decide which of them 
is dominant in the Lorentz force term. We find that the former dominates when n < —1.5, while the latter dominates 
for n > —1.5. However, one cannot determine the relative signs of those two terms unless one also specifies a model 
for the generation of the PMF [sj- Thus, we must decompose the factors into various combinations, i.e. 

S = S[LF] X S[DF], (44) 

where 

S[LF] = < -I, n> -1.5 (II), (45) 



and 



< S[DF] < 1 (i), 
S[DF] = (ii), 
-1 < S[DF] < (iii). 



In the different regimes, S[lf] represents either: (I) the pressure dominated case; (II) the tension dominated case, 
where the magnetic field pressure and tension forces act in the same direction; or (III) the tension dominated case, 
where the magnetic field pressure and tension forces act in the opposite direction. On the other hand, sppj represents 
either: (i) a positive correlation between the matter and PMF distributions; (ii) no correlation; or (iii) a negative 
correlation. Thus, if s < 0, the matter and PMF distributions could be correlated positively (spF] > 0) and the PMF 
pressure dominates in the Lorentz term (for n < —1.5). Another possibility is that the matter and PMF distributions 
negatively correlate (spF] < 0) and the PMF tension dominates in the Lorentz term (for n > —1.5) and the tension 
acts on the density field in the same direction as the magnetic field pressure. Yet another possibility is that the PMF 
tension dominates in the Lorentz term (n > —1.5), but the tension acts on the density field in the opposite direction 
from the pressure force. In these cases the PMF effects act like a gas pressure to oppose the gravitational collapse 
and cause the density perturbations to more slowly evolve. 

On the other hand, if s > 0, the matter and PMF distributions could positively correlate (spF] > 0) and the PMF 
tension dominates in the Lorentz term {n > —1.5) while the tension acts on the density field in the opposite direction 
from the pressure force. Alternatively, the matter and PMF distributions could negatively correlate (s[df] < 0) and 
the PMF pressure dominate in the Lorentz term (n < —1.5). In these cases the Lorentz force from the PMF accelerates 
the gravitational collapse. After decoupling, S does not oscillate and the perturbation evolution is straightforward for 
all of the above cases. 



V. RESULTS AND DISCUSSIONS 



For completeness, in this section we briefiy review the effects of the PMF on the cosmological density field fiuctua- 
tions (see |33| for details). We will then illustrate that relation between erg and the PMF parameters. We will show 
that the constraints on ag from observation give a strong constraint on the PMF parameters. Since parameters of 
the PMF have a strong degeneracy, the existence of such a prior, can be used to effectively constrain the PMF. Also, 
since as is constrained by diverse observational data on linear cosmological scales, we can obtain a reliable prior for 
use in determining likelihood functions for the parameters of the PMF from CMB observations. The PMF effects 
dominate the matter power spectra for wavenumbers k > 0.1 Mpc~^ [STj. This is because the PMF energy density 
fluctuations depend only on the scale factor a and can survive below the Silk damping scale. Therefore, the PMF 
continues to source the fluctuations through the Lorentz force even below the Silk damping scale. In the case of no 
correlation between the PMF and the matter density fluctuations, the matter power spectrum is increased by the 
PMF, independently of whether the PMF pressure or tension dominates. 

Here we note the different effects of the PMF on the power spectrum function P{k) and the matter density fluctuation 
i5. While the total density fluctuation 6 can be smaller or larger depending upon whether the effect of the PMF is 
dominated by its pressure or tension, the power spectrum function P{k) always increases when the PMF does not 
correlate with the primordial density fluctuations. This is because P{k) cx and is not affected by the sign of S. 



A. Effects of PMF parameters on as 



The alternative normalization parameter as is the root-mean-square of the matter density fluctuation in a comoving 
sphere of radius 8h~^ Mpc. It is given by a weighted integral of the matter power spectrum [42]. We can study the 
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physical processes of density field fluctuations on cosmological scales within the linear regime to determine erg. Recently 
CTg has been constrained by observations [H, [s^, [g^, [6l| to be in the range 0.7 < erg < 0.9. From this we can obtain 
strong constraint for the PMF parameters by numerically calculating ug under the influence of PMF effects. 

We expect that the discrepancy between theoretical estimates and observational temperature fluctuations of the 
CMB for higher multipolarity {£ > 1000) is solved by combining a PMF of strength 2.0 nG < \Bx\ < 3.0 nG and the 
SZ effects. In this case, ag derived by such a field strength for the PMF is 0.77 — 0.88. This is consistent with our 
assumed prior in the range erg as 0.7 < erg < 0.9. Since erg is affected by other cosmological parameters, fit, ^^cdm, ns, 
and As, we should consider the degeneracy between the PMF and other cosmological parameters as mentioned above. 
Fortunately, these cosmological parameters are constrained by recent CMB observations on larger scales {£ < 1000) 
[illel,!!!!, while it was shown in our previous work [2ll,[2l|33 that the effect of the PMF mainly affects the CMB 
anisotropics on smaller scales { i > 1000). Hence, we expect that the degeneracy between the PMF parameters and 
the other cosmological parameters is small. For this reason in the present analysis we are justified in fixing the other 
cosmological parameters at their best fit values. 

Figure[l]shows the behavior of the PMF parameters B\ and ne for various constant values of erg as labeled. Since the 
PMF power spectrum depends upon ub from Eqs. (fTSlfTS]) [25l| . for ub > —1.5, the PMF effects on density fluctuations 
for the small scales decrease with lower values for ub- While, for ub < —1.5, P(fc)pMF oc 5|^(^"B+3)/(nB+5)] ^^^^ 
Eqs. p]fT5l)|25t. Therefore, for a spectral index near tt-b = —3.0, the matter power spectrum in the presence of a 
PMF is smaller, and larger amplitudes of B\ are allowed. Since the correlation effects of the PMF for negative and 
positive correlations change at ub — —1.5 (see Section WV} . we divide the discussion below of each correlation into 
two parts based upon whether the spectral index is greater or less than ub = —1.5. 

1. No Correlation 

When there is no correlation between the PMF and the density fluctuations from primary perturbations, the third 
terms in both Eqs. (|42)l and ((43)) vanish and the PMF only acts to increases the total matter power spectrum. In 
this paper, we adopt the constraint that PMF parameters giving erg > 1 are excluded by observations. Panel (a) of 
Figured] shows that a PMF amplitude of B\ ^1 nG is excluded when ub > —0.9. Furthermore, PMF amplitudes 
of B\ ^0.11 nG are excluded when tt-b > 0.2. The magnetic field strength in galaxy clusters is ~ 1/iG. Therefore, 
if isotropic collapse is the only process which amplifies the magnetic field strength, the lower limit to the PMF is 
^ InG at z ^ 0. Hence, we can obtain a strong constraint on this PMF evolution model for a PMF spectral index 
in the range ub < —0.9. 

2. Negative Correlation 

For the case of negative correlations, the pressure of the PMF dominates for ub < —1.5, and the PMF causes an 
increase in the density fluctuations. For ub > —1.5, however, the tension of the PMF dominates, and the PMF causes 
a decrease in the density fluctuations. These behaviors can be traced to the third terms in each Eqs. (|l^ and (I43p . 
Using the allowed range of PMF parameters as mentioned above, a PMF of B\ ^1 nG is excluded for ub > —0.81, 
and a PMF of Bx >0.11 nC is excluded for ub > 0.26. 

3. Positive Correlation 

When there is a positive correlation, the pressure of the PMF dominates for ub < —1.5, and the PMF leads to a 
decrease in the density fluctuations. On the other hand, the tension of the PMF dominates for nB > —1.5. In this 
case the PMF causes an increase in the density fluctuations. We can attribute these behaviors again to the third 
terms in Eqs. and Using the allowed range of PMF parameters as noted above, a PMF of B\ ^1 nG is 

excluded for ub > -0.94, and a PMF of B^ >0.11 nG is excluded for ub > 0.13 

For both negative correlations with ^b > —1.5 and positive correlations with nB < —1.5, the PMF decreases 
the total matter power spectrum until the strength of the PMF effect is comparable to the primary matter power 
spectrum. Beyond this point, the density fluctuations from the PMF exceed the density fluctuations from the primary 
power spectrum and the PMF effect dominates the total matter power spectrum. The amplitude of the total matter 
power spectra when a PMF is present is greater than the total matter power spectrum without the presence of a 
PMF. 
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In these cases, we obtain different constraints on the strength of the PMF for each value of Ub from the other 
models. We assume that ranges of the PMF parameters giving erg < 0.6 are excluded by observations. Panel (c-1) of 
Figure H] shows the excluded range of the PMF parameters within the contour defined by erg = 0.6. 

Our result is consistent with previous constraints on PMF parameters (22l . [56j , and our new more precise constraints 
obtained from the matter power spectrum are independent of previous methods to constrain the PMF parameters. 
Hence, we can now constrain more precisely the physical processes by which a PMF affects the evolution of structure 
on cosmological scales. 



VI. SUMMARY 



A primordial magnetic field (PMF) affects the evolution of density field fluctuations in the early universe. Therefore, 
we can constrain PMF parameters, e.g. the PMF amplitude Bx and power spectral index riB, by comparing a 
theoretical calculation of the density field fluctuations affected by a PMF and observational data, e.g. a number 
density fluctuation of galaxies as indicated by the cg parameter. We have illustrated the relation between erg and the 
PMF parameters. We have shown that the observed range of cg gives a strong limit on the PMF parameters given 
by riB < 0.9 (no correlation case), < —0.86 (negative correlation case), and < 0.94 (positive correlation case) for 
Bx > InG. We have also shown how the PMF parameters which were constrained by previous methods are affected 
by our new constraint from erg. Since density field fluctuations are the origin of the LSS on cosmological scales, we can 
use this as a prior to study the physical processes of the PMF and to place better constraints on the PMF parameters. 
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